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Abstract. 

The objective of this introduction to Colombeau algebras of generalized-functions (in 
which distributions can be freely multiplied) is to explain in elementary terms the essential 
concepts necessary for their application to basic non-linear problems in classical physics, 
ly-j I Examples are given in hydrodynamics and electrodynamics. The problem of the self- 

J> ' energy of a point electric charge is worked out in detail: The Coulomb potential and 

. field are defined as Colombeau generalized-functions, and integrals of nonlinear expressions 

I corresponding to products of distributions (such as the square of the Coulomb field and the 

square of the delta-function) are calculated. 

Finally, the methods introduced in Eur. J. Phys. 28 (2007) 267, 1021, and 1241, to deal 
' with point-like singularities in classical electrodynamics are confirmed. 

> 

a 
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1. Introduction 



The theory of distributions, invented by Laurent Schwartz nearly 60 years ago, provided a 
simple and rigorous calculus unifying a great variety of previously ill-defined mathematical 
^ ' techniques used in physics and engineering. Improper functions such as Heaviside's step 
I function H (x) and its derivative, Dirac ' s delta function 6{x), were replaced by linear functionals 
called 'distributions,' which like C°° functions can be differentiate any number of times III 13- 
Distributions, however, cannot in general be multiplied. In mathematical language the set 
V of all distributions is a vector space rather than an algebra. But products of distributions arise 
naturally in many areas of science and engineering, most prominently in electrodynamics and 
particle physics as self-interaction terms of the type and and in hydrodynamics 

as products like H(x)(5(x) of shock waves and their derivatives, as well as in mathematics as 
solutions of partial differential equations. Many proposals have therefore been made to define 
an algebra of generalized functions Q such that V c G. 

The difficulty of this task is not only mathematical (Schwartz's theory of distribution is 
highly abstract) but also conceptual since many generalizations are possible. In this sense the 
algebras defined by Jean-Francois Colombeau (SHHHl have essentially optimal properties for 
a wide range of applications Il6l|71[8l|9|, which combined with the fact that they provide what 
is possibly the most simple and natural generalization of the space of distributions, have made 
them very popular (Ref. [lOJ and numerous references therein). 
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Trying to make these developments accessible to as wide an audience as possible is also 
a difficult task: This paper is therefore a rather specialized article, at the crossroad of modem 
mathematics and physics, which would only be accessible to advanced undergraduates and 
graduate students (and their teachers) in mathematical and theoretical physics. 

The ambition of this paper is to give a straightforward introduction to Colombeau 
generalized functions and to their applications by emphasizing the underlying concepts rather 
than the mathematical details of the theory, in order to make it accessible to most physicists 
and engineers. The idea is that working with Colombeau functions is like working with real 
numbers: There is no need to know how transcendental numbers are embedded in the set of 
real numbers to use them. 

In particular, we do not introduce Colombeau's algebras by starting with their abstract 
definition as quotient spaces Q = £m/J^, where £m is a space of 'moderate' (or 'multipliable') 
functions, and A/" C an ideal of 'negligible' functions such that £m/J^ becomes a superset 
of v. We prefer to emphasize the fundamental idea behind their construction, which is that the 
smooth functions (i.e., the C°° functions — which are indefinitely continuously differentiable) 
must be a faithful subalgebra of Q. Moreover, to be self-consistent, we begin by recalling in 
Sec.[2]some basic definitions, and in Sec.[3]the fundamentals of Schwartz distributions. In the 
same spirit we summarize in Sec. HI that is before defining Colombeau's algebra in Sec.Hl the 
reasons why the product of distributions is inconsistent, and how regularization can be used 
to circumvent this problem. 

We then consider two sets of applications. First, in Sec.|71 to hydrodynamics, the second 
major topic after quantum field theory to which Colombeau originally applied his theory 
llTl, which enables to illustrate its power in numerical engineering and as a general method 
for solving nonlinear partial differential equations. Second, in Sees. [8] to \Wi to classical 
electrodynamics, which enables to show how the introduction of point charges using the 
Colombeau formalism leads to the possibility of calculating divergent quantities such as 
the self-energy — which are quadratic in the fields, and therefore meaningless in Schwartz 
distribution theory. Finally, in Sec. [HI we use the formalism of generalized functions to 
confirm the simplified methods introduced in References (TV\ and [ | 12| to deal with point-like 
singularities in linear problems of classical electrodynamics. 

2. Definitions and notations 

This section is put here rather than in appendix because it recalls definitions and notions which 
give an opportunity to anticipate important concepts that will be developed in the text, and 
because the subject of generalized functions makes it essential to be particularly careful with 
notations, and with the belonging of all objects to the proper sets. 

fl An open set of M". In this paper f2 is any interval ]a, b[c M such that a < < 6, unless 
otherwise specified. 

No The set {0, N}, i.e., the natural numbers and zero. 
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C°° The algebra of continuous functions on Vt which are smooth, i.e., that are indefinitely 
continuously differentiable, and which have compact support. 

The space of continuous functions on which are m-times continuously differentiable 
and which have compact support. 
C The space C = C° of continuous functions on Vt with compact support. Such functions 

may have points at which the left and right derivatives are different. 
Cp The space of piecewise continuous functions on Vt with compact support, i.e., continuous 
except on a discrete set, on which they have left and right limits, and where differentiation 
leads to 5-functions. 

V The space of functions C^iVt) equipped with an inductive limit topology suitable to 
define Schwartz distributions. The functions in this space are used to construct the 
weakly converging sequencies defining distributions, as well as to provide so-called 'test 
functions,' denoted in this paper by T, on which the distributions are evaluated. 

V The space of Schwartz distributions, i.e., the set of all linear continuous functions (linear 
forms) on V{Vl), that is, the topological dual of V{VL). 

Q The algebra of Colombeau generalized functions, or ^-functions, on Vl. 

The Colombeau algebra Q, as well as three function spaces: £,£^, and M, will be 
defined in Sec. [H 

In the sequel we will generally keep n implicit, and tacitly assume that all functions are 
zero outside of f2. Finally, we will use the term 'sequence' for expressions such as lim^^o Fe 
even though mathematicians reserve this term to mappings n i— > F„, so that e ^ corresponds 

to 1/n ^ as n ^ oo. 

3. Schwartz distributions 

The space V of Schwartz distributions contains, besides the ordinary (i.e., C, and C°°) 
functions, generalized functions corresponding to discontinuous functions and unbounded 
functions. While these functions cannot be differentiated in the classical sense, they can be 
indefinitely differentiated in the sense of distributions. That is, if D E V is any distribution, 
its derivatives in the 'distributional sense' are such that, VT G V and D" = d'^/dx'^. 



It is not possible to represented non-trivial distributions, such as Dirac's 5-function, by 
simple algebraic formulas or even by ordinary limiting processes. They can however be 
represented by sequences of smooth functions 



for which ordinary pointwise convergence is not required. Instead, what is required is 'weak 
convergence' for the scalar product of D{x) with any test functions T{x) G V, i.e., the 
existence of the limit 




(3.1) 



D{x) := limD^{x), where D^{x) G V, 



(3.2) 



VT G V, 




(3.3) 
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The meaning of operating 'in the sense of distributions' is then that all operations on 
distributions are actually performed on D^(x), while D(x) can be seen as a convenient symbol 
to designate a given distribution. 

Eq. (13.31) shows that distributions can be interpreted as linear functionals D(T) = {{D\T)) 
defined by their effect on test functions. Moreover, since many different sequencies may 
converge weakly to the same limit, each distribution corresponds to an equivalence classes of 
such sequencies, which all together form the Schwartz distribution space V: 

Definition 1 Two distributions D and E G V, of respective representatives and 
E^, are said to be equal (or equivalent), and one write D = E, iff 

VT G V, lim j {d,{x) - E,{x)^T{x) dx = 0. (3.4) 

For instance, Dirac's 5-f unction is defined by the property 

{{S\T)) = [ 5{x)T{x) dx = T(0), (3.5) 
Jn 

so that all sequences which have this property form an equivalence class corresponding to 
Dirac's 5-function distribution, conventionally denoted by the symbol '(5(s).' Two examples 
of such sequences are 

1 1 
6,{x) = — — and 6,{x) = ^^exp{ — -). (3.6) 

More generally, any normalizable C°° function p(y) with compact support can be used to 
define (5-sequencies, i.e., 

WpeC^, [p{y)dy = l 6,{x):=-p{-). (3.7) 

The requirement that both and T are functions and the definition (13.11) enable 

to derive at once a number of useful properties. For instance, the equations 

{{x6{x)\T{x))) =0, and {{-6{x)\T{x))) = -T'{0), (3.8) 

X 

which are often symbolically written 'x5{x) = 0,' and 'x^^5{x) = — 5'(x)' or 'xS'{x) = 
—6{xy are the fundamental formulas of calculus with distributions Q 

In summary distributions are not functions in the usual sense but equivalence classes 
of weakly convergent sequencies of smooth functions. All operations on distributions are 
therefore made on these sequencies, which are thus added, differentiated, etc., according to 
the operation in question. It is remarkable that distributions enjoy essentially all properties 
of C°° functions, including multiplication by a C°° function, with a few exceptions such as 
the impossibility to multiply two distributions in the general case, as was demonstrated by 
Schwartz in his impossibility theorem of 1954 [|71 p. 8], IfTOl p. 6]. 

Moreover, Schwartz distributions have a very precise relation to continuous functions, 
which can be spelled in the form of the theorem: 

^ As will be seen, the Colombeau formalism provides a non-ambiguous notation for these formulas, i.e., 

x6{x) X 0, and x^^S{x) x —S'{x) or xS'{x) x —S{x). 
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Theorem 1 (Schwartz local structure theorem) Any distribution is locally a partial 
derivative of a continuous function [|7l p. 6]. 

Differentiation induces therefore the following remarkable cascade of relationships: 
continuously differentiable functions continuous functions distributions. This gives 
a unique position to Schwartz distributions because they constitute the smallest space in which 
all continuous functions can be differentiated any number of times. For this reason it is best 
to reserve the term 'distribution' to them, and to use the expression 'generalized function' for 
any of their generalizations. 

4. Multiplication and regularization of distributions 

There are two kinds of problems with the multiplication of distributions: (i) The product 
of two distributions is, in general, not defined. For example, the square of Dirac's 5- 
function is not a weakly converging sequence, as can easily be verified by squaring either 
of the sequencies in Eq. (13.61) and trying to evaluate them in a scalar product with any test 
function, (ii) Differentiation is inconsistent with multiplication because the Leibniz rule, or 
even associativity, can fail under various circumstances. For example, while Dirac's 5-function 
is related to Heaviside's step function through differentiation as 6{x) = H'(a;), the algebraic 
identity H^(x) = H(x) leads to inconsistencies. Indeed, 



Over the years many methods for solving these problems have been proposed. One of the 
simplest and most effective is 'regularization,' which consists of modifying the functions to 
be multiplied or differentiated in such a way that they become more regular (i.e., continuous, 
differentiable, finite, etc.). All operations are then done with the regularized functions until 
the end of the calculation, and the final result is obtained by the inverse process which returns 
the function from its regularization. 

A particularly convenient regularization technique is based on the convolution product. 
For instance, if f{x) is any function on R, its regularization f^{x) is 



Here Pe{x) is a smoothing kernel (also called regularizer or 'moUifier') which in its simplest 
form is a 5-sequence as defined in Eq. (|3.7I) . and e g]0, 1[ is the (so called) regularization 
parameter. Consequently, in the limit e 0, the moUifier becomes equal to the (5-function, 
which by Eq. (13.51) acts as the unit element in the convolution product, i.e., f * S = f. Thus, 
when e 7^ the regularization is such that / is 'mollified' by the convolution, while / can be 
retrieved by taking the limit e ^ 0. 

The power of convolution as a regularization technique stems from the theorem: 

Theorem 2 The convolution {D * p) (x) of a distribution D eV by a function p eV 
is a C°° function in the variable a; [[B P- 465]. 



= H ^ 2m = 5 2Y(^5 = H5 ^ 2H5 = H5 ! 



(4.1) 




(4.2) 



Colombeau generalized-functions in classical electrodynamics 



6 



Regularized functions f * Pe and distributions D * can therefore be freely multiplied and 
differentiated. Moreover, the mollified sequence = D * provides a representative 
sequence of the type (|3.2I) of any distribution D E V. 



5. Colombeau generalized functions 

A Colombeau algebra Q is an associative differential algebra in which multiplication, 
differentiation, and integration are similar to those of C°° functions. Colombeau and others 
have introduced a number of variants of Q but all 'Colombeau algebras' have in common one 
essential feature: The C°° functions are a faithful differential subalgebra of Q, a feature that 
Colombeau discovered to be essential to overcome Schwartz's multiplication-impossibility 
theorem. 

With hindsight it is easy to understand why: If we suppose that Q is an algebra containing 
the distributions and such that all its elements can be freely multiplied and differentiated just 
like C°° functions (i.e., in a way respecting commutativity, associativity, and the Leibniz rule), 
then C°° must be a subalgebra of Q because C°° C V. Thus, to define Q, it suffices to start 
from a differential algebra £ containing the distributions, and then to define ^ as a subalgebra 
of £ such that the embedding of the C°° functions in Q is an identity. In formulas: if [g] E Q 
represents an object g embedded in Q, then for all / E we want that [/] = /, whereas for 
any other function or distribution D E V we may have [D] ^ D. 

This simple observation gives a powerful hint for an elementary construction of Q because 
by Definition |3] there is a one to one correspondence between any arbitrary distribution D{x) 
and a class of weakly convergent sequence of C°° functions D^(x), and by Theorem |4] any 
representative of that class can be written as a convolution of the form (|4.2I) . Thus, the starting 
point is to consider for S the set of mollified sequencie^ 

S:=\f,:{r],x)^Ur],x)}, (5.1) 



which are C°° functions in the variable x for any given Colombeau mollifier r], and depend on 
the parameter e g]0, 1[ through the scaled moUifieiE 

?7e(a;) := -?7(^— j , normalized as j dyr]{y) = l. (5.2) 

The distributions f E V are then embedded in S as the convolutio 

1 fy — X 



fe{x) := * f{x) = j dy -^(^-^) f{y) 

dz ri{z) f{x + ez), (5.3) 



where, in order to define Q C S, the Colombeau mollifier i] may need to have specific 
properties in addition to the normalization (|5.2I) . 

^ In the literature the notation {fc)c is often used here instead of f^. It emphasizes that is an element of £ 
rather than just a representative sequence or a regularization. 

^ In this and the next sections we take 17 = M so that all integrals are ■ The generaUzation to fi = M'^ is 
immediate, e.g., the scaled mollifier is rj{x/e)/e^. 

^ This definition due to Colombeau differs by a sign from the usual definition ( I4.2l i of regularization. 
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To find these additional properties we have to study the embeddings of functions. We 
therefore calculate (15.31) for / G C°°, which enables to apply Taylor's theorem with remainder 
to obtain at once 



Ux)=f{x) J dzr]iz) + ... (5.4) 

+ / dz z^'viz) + ... (5.5) 

n\ J 

+ -^^^ j dz z'^^'r^iz) f^^+'\x + ^ez), (5.6) 

where /("^(s) is the ra-th derivative of f{x), and i!) g]0, 1[. Since / G C°° and r] has compact 
support, the integral in (15.61) is bounded so that the remainder is of order 0(e''^^) at any fixed 
point X. 

Then, if following Colombeau the moUifier 7] is chosen in the set 

and /d..",W^O. V„^l.....,eN}. (5.7) 

all the terms in (15.51) with n E [1, g] are zero and we are left with 

V/gC~, /,(x)=/(x) + 0(e'^+^). (5.8) 

Therefore, provided r/ is a Colombeau moUifier and q can take any value in N, it is possible to 
make the difference feix) — f{x) as small as we please even if e g]0, 1 [ is kept finite. 

In terms of the embeddings (15.31) the condition [/] = / insuring that the smooth functions 
are identically embedded in Q is that [/e](a;) = [f]{x) = f{x) for all / G C°°. Thus, 
comparing with (|5.8I) . we are led to consider the set M of the so-called negligible functions, 
which correspond to the differences between the functions and their embeddings in S, i.e., 

V/gC°°, VgGN, f,{x)-f{x) = 0{e'') eAf. (5.9) 

To define Q we need a prescription such that the differences (15.91) can be neglected, i.e., 
equated to zero in Q. Moreover, for Q to be an algebra, that prescription must be stable under 
multiplication. That means that all elements g E Q have to have the property that any of their 
representatives ge ^ £ multiplied by a negligible function are negligible. In mathematical 
language, J\f has to be an ideal of the subset {g^} = £m C £ of all representatives of all 
elements of Q. Or, in simple language, the negligible functions have to behave as the 'function 
zero' when multiplying any function of £m. It is however very simple to characterize this 
subset: Following Colombeau we call the elements of £m moderate ( or multipliable ) functions, 
and we define 

\/g,E£M: 3iVGNo, such that g,{x) = 0{e-^). (5.10) 

Indeed, as q in (15.81) is as large as we please, and N in (15.101) a fixed integer, the product of a 
negligible function by a moderate one will always be a negligible function: M is an ideal of 
£m. Moreover, the product of two moderate functions is still moderate: They are multipliable. 



For example, the Colombeau embeddings (15.31) of the 5 and Heaviside functions are 

I / x\ r^^ 

5^(a;) = -r/( — ), and B.,{x) = dzri{~z), (5.11) 



e 



— oo 
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which are moderate functions with = 1 and 0, respectively. More generally, it can easily 
be proved using Schwartz's local structure theorem that: 

Theorem 3 (Colombeau local structure theorem) Any distribution is locally a 
moderate (i.e., multipliable) generalized function [[31 p. 61], [iTOl p. 62]. 

Therefore, Af C C (C)^ C {T>')^ C £m- It is also a matter of elementary calculation to 

verify that and J\f are algebras for the usual pointwise operations in £. Moreover, £m is a 
differential algebra, and it is not difficult to show that £m is the largest differential subalgebra 
(i.e., stable under partial differentiation) of £ in which J\f is a differential ideal. 

The fact that J\f is an ideal of £m is the key to defining Q. Indeed, if we conventionally 
write J\f for any negligible function, then 

V(7„ h, G £m, (g. + AT) • (h, +M)=g,-K + U. (5.12) 

Thus, it suffices to define the elements of Q as the elements of £u modulo M, i.e., to define 
the Colombeau algebra as the quotient 

Q ■■= % (5.13) 

That is, an element 5? G ^ is an equivalence class [g] = [g^ + A/"] of an element g^ G £m, which 
is called a representative of the generalized function g. If '0' denotes multiplication in Q, the 
product g Qhis defined as the class of g^ ■ where g^ and are (arbitrary) representatives of 
g and h; similarly Dg is the class of Dg^ if D is any partial differentiation operator. Therefore, 
when working in Q, all algebraic and differential operations (as well as composition of 
functions, etc.) are performed component-wise at the level of the representatives g^- 

Q is an associative and commutative differential algebra because both £m and J\f are such. 
The two main ingredients which led to its definition are the primacy given to functions, 
and the use of Colombeau moUifiers for the embeddings. In fact, Colombeau proved that the 



set (15.71) is not empty and provided a recursive algorithm for constructing the corresponding 
moUifiers for all g G N. He also showed [|7l p. 169] that the Fourier transformation provides 
a simple characterization of the moUifiers. But, in the present paper as in most applications 
of the Colombeau algebras, the explicit knowledge of the form of the Colombeau moUifiers is 



not necessary: It is sufficient to know their defining properties (15.71) . 

For example, let us verify that 6e{x) given by (|5.11l) has indeed the sifting property 
expected for Dirac's (5-function. Starting from (13.11) and employing Taylor's theorem we can 
write 



{{S,\T)) = f 6,{x)T{x) dx = f dx-7](--)T{x)= f dz r]{-z)T{ez) 

J J 6X6/ J 



= dz vi-z) (t(0) + ezT'iO) + ^T"(0) + ...), (5.14) 

Then, in Schwartz theory, we take the limit (13.31) . i.e., 

lim ((5,|r)) = r(0) + O(6), (5.15) 

which is the expected result thanks to the normalization (15.21) . However, in Colombeau theory, 
there is no need to take a limit to get the sifting property because in the development (15.141 ) 
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the conditions (15.71) imply that all terms in with l<n<g+lare identically zero. Thus 

((4|T))=T(0) + O(e''+i), VgGN, (5.16) 

where the remainder is an element of J\f so that in Q the sifting property of 6^ is an equality 
rather than a limit. It is this kind of qualitative difference between the Schwartz and Colombeau 
theories which makes it possible in Q to go beyond distribution theory. 



6. Interpretation and multiplication of distributions 

To construct the Colombeau algebra we have been led to embed the distributions as the 



representative sequences 7e G £^ defined by (|5.3I) where r] is a Colombeau moUifier (15.71) . We 
can therefore recover any distribution 7 by means of (13.31) . i.e., as the equivalence class 



7(T) := lim j dx -f,{x) T{x), Vr(x) G V, (6.1) 

where 7^ can be any representative of the class [7] = [7^ + A/"] because negligible elements 
are zero in the limit e — > 0. 

Of course, as we work in Q and its elements get algebraically combined with other 
elements, there can be generalized functions [g^:] different from the class [7^] of an embedded 
distribution which nevertheless correspond to the same distribution 7. This leads to the concept 
of association, which is defined as follows^ 

Definition 2 Two generalized functions g and h E Q, of respective representatives 
g^ and h^, are said to be associated, and one write g ^ h, iff 

lim j dx (g,{x) - h,{x)^T{x) = 0, VT(x) G V. (6.2) 

Thus, if (7 is a generalized function and 7 a distribution, the relation (7x7 implies that g admits 
7 as 'associated distribution,' and 7 is called the 'distributional shadow' (or 'distributional 
projection') of g. 

Objects (functions, numbers, etc.) which are equivalent to zero in Q, i.e., equal to 
O(e^), Vg G N, are called 'zero.' On the other hand, objects associated to zero in Q, that is 
which tend to zero as e ^ 0, are called 'infinitesimals.' Definition (|6.2I) therefore means that 
two different generalized functions associated to the same distribution differ by an infinitesimal. 

The space of distributions is not a subalgebra of Q. Thus we do not normally expect that 
the product of two distributions in Q will be associated to a third distribution: In general their 
product will be a genuine generalized function. 

For example, the square of Dirac's 5-function, Eq. (15.1 II) . which corresponds to 
(5^)e(x) = {5^y{x) = e'^'^rf{—x/e), has no associated distribution. Indeed, making a Taylor 
development as in (|5.14l) . 

lim((5^|T)) = lim / dx—rfi )T(x) = lim- / dz rfi—z^Tiez) 

^ In the literature the symbol « is generally used for association. We prefer to use x because association is 
not some kind of an 'approximate' relationship, but rather the precise statement that a generalized function 
corresponds to a distribution. 
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= Jim 1 dz r^\-z) (t(0) + ezT'(O) + ^T"(0) + ...) 

= lim^^ /" dzr]\~z)+T'{0) [ dz zrf{-z) = oo. (63) 

But, referring to (15.101) . is a moderate function with N = 2. The square of 5{x) makes 

therefore sense in ^ as a 'generalized function' with representative = rf{—x/e)/e^. 

Moreover, its point- value at zero, ?7^(0)/e^, can be considered as a 'generalized number.' 

On the other hand, we have in Q elements like the ra-th power of Heaviside's function, 
Eq. (15.1 II) . which has an associated distribution but is such that [H"] (x) 7^ [H] {%) in Q, whereas 
H"(a;) = H(x) as a distribution in V . Similarly, we have [x] 7^ in ^, whereas 

x5{x) = in V . In both cases everything is consistent: Using (|6.2I) one easily verifies that 
indeed [H'^](x) x [H](a;) and [x] [5]{x) x 0. 

These differences between products in Q and in V stem from the fact that distributions 
embedded and multiplied in Q carry along with them infinitesimal information on their 
'microscopic structure.' That information is necessary in order that the products and their 
derivatives are well defined in Q, and is lost when the factors are identified with their 
distributional projection in V. For example, since [H"](x) 7^ [H](a:;) the inconsistencies 
displayed in (|4.1I) do not arise in Q. Nevertheless, if at some point of a calculation it is 
desirable to look at the intermediate results from the point of view of distribution theory, one 
can always use the concept of association to retrieve their distributional content. In fact, this 
is facilitated by a few simple formulas which easily derive from the definition (16.21) . For 
instance, 

V/i,V/2eC ^ [/i] [/2] X [/i • /2], (6.4) 

V/GC°°,V7eP' ^ [f]&[l]-[f -li (6.5) 

V71, V72 e ^ [71] [72] ^ [71 ■ 72], (6.6) 

^91,^92 eg, 91-92 D"(7ixD-(72. (6.7) 

For example, applying the last equation to [H^](a;) x [H](x) one proves the often used 
distributional identity 2 [5] (2;)[H] (a;) x [6]{x). 

In summary, one calculates in Q as in C°° by operating on the representatives 9e ^ £ with 
the usual operations {+, — , x , d/dx}. The distributional aspects, if required, can be retrieved 
at all stages by means of association. However, as will be seen in the applications, it is possible 
in many cases to set aside the concept of distributions and to replace it by the more general 
and flexible one of ^-functions. 



7. Applications to hydrodynamics 

Distributions and their applications are common place in most areas of physics and engineering. 
It is well known how to evaluate a (5-function, how to calculate with piecewise continuous 
functions, etc. In such applications there is no essential difference between Schwartz 
distributions and Colombeau functions. This is because most characteristics distinguishing 
(5-sequencies like for example those of Eq. (|3.6I) have no effect in such applications since 
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the only things that matter when they are evaluated on test functions is their sifting property 
and their integral, which is normalized to one. These characteristics are however relevant 
when distributions are multiplied and evaluated at the same x in expressions such as 5'^{x), 
H(x)5(x), Hi(x)H2(x), 5i{x)52{x), Y{i{x)52{x), etc., where the indices 1 and 2 refer to 
different 5-functions such that Y{[ = 5i. 

For instance, products of the type Hi(a;)H2(a;) and }li{x)52{x) arise in the study of 
the propagation and interaction of shock waves, such as those occurring in strong collisions 
between projectiles and armor, a subject that has been extensively studies by Colombeau and 
his collaborators (Ref. |I71 and references therein). Shock waves induce sudden and large 
variations of physical quantities, for example the density, on a distance comparable to only 
a few times the average distance between molecules. An Heaviside step function would 
therefore appear to be an excellent approximation of that behavior, which is indeed often the 
case. However, in strong shocks during which a phase transition (e.g., from elastic to plastic) 
occurs this approximation is insufficient. For instance, a typical combination of distributions 
arising in such a case are products like 

Hi(x)H2(x) X a5(x), (7.1) 

where 5 is H'^ or Hg. It then turns out that whereas a = 1/2 in the simple case Hi = H2, 
measurements show that a can be anything between 0.05 and 0.95, see \Jj p. 43-48], implying 
that the 'microscopic profiles' of Hi and H2 are very different at the jump, i.e., 5i 7^ 62- 

A lesson from this application is that by assuming that physically relevant distributions 
such as H and 5 are elements of Q one gets a picture that is much closer to reality than if they 
are restricted to V. In fact, this lesson applies not just to numerical modeling and applied 
physics but also to mathematics and theoretical physics. 

Consider for example one of the simplest nonlinear partial differential equations, the 
inviscid Burger's equation of hydrodynamics, which in Q can be written in two ways 

Ut + uux = 0, or ut + uux X 0. (7.2) 

Both equations have a traveling wave solutions of the type u{x, t) = {u2 — ui)Y{{x — ct) + ui, 
which using the identity 2HH' = H' yields a jump velocity of c = {ui + U2)/2. But, if 
multiplied by u, the first equation has an additional solution which turns out to be inconsistent 
with the first one, so that it has in fact no solutions. This is however not the case with the second 
equation because whereas multiplication is compatible with equality in Q, it is not compatible 
with association. Therefore, the distinction between = and x automatically insures that the 
physically correct solution is selected, a distinction that can be made in analytical as well as 
in numerical calculations by using a suitable algorithm 

8. Point charges in classical electrodynamics 

It is well known that classical electrodynamics is basically a continuum theory. Nevertheless, 
point charges are very often considered — even in the most elementary introductory lectures. 
Whenever a contradiction arises the difficulty is then set aside: Infinite quantities are discarded 
or renormalized, self-interactions terms in Lagrangians are ignored, etc. 
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It is however possible to consistently introduce point charges through distribution theory, 
and to deal with them successfully, at least as long as the electromagnetic fields, currents, 
and charge densities interpreted as distributions are not multiplied^ Since distributions are a 
subspace of Q nothing fundamentally new has to be invented, and the electromagnetic field 
distributions defined in the context of Schwartz distribution theory can just as well be used 
in the ^-setting. This provides a 'standard methodology' in which these distributions are 
simply embedded in Q by means of a Colombeau regularization. But, since the Colombeau 
regularization has specific properties that standard regularizations do not have, it is possible to 
develop a more powerful and convenient methodology which takes these properties explicitly 
into account. This methodology, based on a suitable generalized function T, has already been 
presented in References {\V\ and [fT2ll . albeit with only a simplified and intuitive justification 
of its validity. In the following we are going to give a rigorous mathematical foundation to it, 
starting with a presentation of the standard methodology, and then moving in successive steps 
to the improved one. 

The basic idea in distribution theory is to replace the classical Coulomb potential e/r of 

a point charge by the weak limit of the sequence of distributions [|T3l p. 144], [[21 p. 51], 

g 

0(r) := -Ha(r), where Ha(r) := hmH(r — a), (8.1) 

r ~ ~ fi^o 

where e is the electric charge of an electron at rest at the origin of a polar coordinate system, 
and r = \f\ the modulus of the radius vector^ Consistent with Schwartz's local structure 
theorem, 0(r) is the derivative of e lima^o log('"/'^)H(r — a), a C° function Vr > 0. The 
cut-off a > insures that 0(r) is a well defined piecewise continuous function for all r > 0, 
whereas the classical Coulomb potential e/r is defined only for r > 0. It is then readily 
verified, using (|3.3i) . that (18.11) is a distribution. Indeed, since Vt = M."^, 

VT G V, [[[ d?r (f){r)T{r) = Ane [ dr rT{r) G M, (8.2) 

JJJm.3 Jq 

because T E V has compact support so that the integral is bounded. 

The next step, in order to be able to calculate the field E = — V0 and the charge density 
Anp = V • -E, is to represent by a mollified sequence. Taking for it the Colombeau form 



(15.31) the Coulomb potential is at once embedded in Q as the representative sequence 

e r°° n(z) 

0,(r) = (-HJ (r) = elim / dz Vr > 0. (8.3) 

r - a^o J r + ez 

One then trivialy verifies that the distributional Coulomb potential (18.11 ) can be recovered by 



letting e — » as in (16.21) . i.e., 

0,(r)x-H,(r):=(/.(r), (8.4) 
r ~ 

which reverts to the classical Coulomb potential e/r as a —> 0. 

^ It is also possible to deal in specific cases with problems that are non-linear in the fields. But this requires ad 
hoc prescriptions, such as defining the energy-momentum tensor, which is quadratic in the fields, directly as a 
distribution. See, e.g., MfTSl . 

^ We write a rather than a in subscript to emphasize that a is not an index but the parameter of a limiting 
sequence. 
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Calculating the embedded Coulomb field is now straightforward because the embedded 
potential (eHa/r)^ is C°° in the variable r. It becomes 

E,(f ) = -V0,(r) = elimf T dz -^^^ - -r,(^))u, (8.5) 

where n = Vr is the unit vector in the direction of r. Introducing the notation 

6a(r) := lim6(r — a), so that lim — rji ) = (-6a) (r), (8.6) 

~ a^o a^o ea e a ~ 

this electric field can be written in the more convenient form 

E,{r) = e((^Hj,(r) - (^y,(r))£ (8.7) 

By an appeal to test functions we easily verify that the field is a distribution, and that the 
(5-function in (18.71) gives a nul contribution when evaluated on a test function. Thus 

i,(f)x^H,(r)M:=E(f), (8.8) 

where E(f) is the distributional Coulomb field which in the limit a ^ yields the classical 
Coulomb field ef/r^. Therefore, the distribution E(f) associated to the ^-function E^, i.e., 
its 'shadow' obtained by projecting it on V, does not contain the (5-function contribution on 
the right of (fO). 

To get the Coulomb charge density we have to calculate the divergence of (18.51) . In 
standard distribution theory one would then ignore the term on the right because it corresponds 
to a (5-function which, as we have just seen, gives no contribution when evaluated on a test 
function. However, in Q, this term cannot be ignored if we subsequently calculate quantities 
in which E^ is a factor in a product. Calculating is therefore somewhat laborious, but still 
elementary. It yields, using V ■ m = 2/r, 

i.p.(r) = V . E.(f) = elto( - J^d. - -„( — ) 

+ - 2 r dz'^, + ^,'(— ) ). (8.9) 

This expression can be rewritten in the less cumbersome form 

4vrp.(r) = e(^(lHJ,(r)-2(lHJ,(r) 

+ {-,^a)Xr) - -{-^a)Xr) - ), (8-10) 

r a a - / 

where we have put the two Heaviside-function terms on the first line and the three Dirac- 
function ones on the second. This charge density is of course much more complicated than 
the single three-dimensional (5-function which is associated to it in the standard distributional 
formalism. But it is the correct result, and a typical example of how quickly calculations 
become complicated when the infinitesimal details of distributions are fully taken into account. 
To calculate the distributional shadow associated to (18.101) we remark that the two Heaviside- 
function terms on the first line cancel each other when e ^ 0, and that, in this limit, the 
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representatives of the 5-functions on the second line become genuine (5-functions. Thus, 
evaluated on a test function T, the first two 5-terms give 

POO 12 1 

/ drr^i— )5{r-a)T{r) = -Tia), (8.11) 

Jq ra 

whereas the S'-term gives, using integration by parts, 

dr )5'{r — a)T{r) = — dr ( T(r)) 5{r — a) 

a Jq a 

= 2T{a) + aT'{a). (8.12) 

Thus, adding (18.1 II) and (|8.12l) . and passing to the limit a ^ 0, we get the test-function- 
evaluation T(0) so that 

P.Wx^5,(r):=p(r), (8.13) 

which yields the classical point-charge density e5^{r) = e5{r)/ATir'^ as a ^ 0. 

9. Point charges in T -methodology 

As we have seen in the previous section, calculations with distributions embedded in Q can 
easily become tedious. However, this is mainly because the representatives used for 0, E, and 
p are mollified embeddings of sequencies of distributions, such as (|8.3I) for </>, in which the 



full properties (15.71 ) of the Colombeau moUifiers have not yet been taken into account. Indeed, 
instead of (18.31) . the Coulomb potential can actually be written in terms of a ^-function T so 
that it reads Dm 

e 1 f for r = 0, 

0(r) = -T(r), where -T(r) = < (9.1) 
r r 1^ r ^ for r > 0, 

which has the advantage that most calculations can be made as if 1/r and T(r) were ordinary 
functions. 

In fact, a representative of the ^-function T is provided by the embedding of the 
distribution, i.e., 

noo 

T(r) := hmUm(Ha) (r) = limlim / dzr]{z), (9.2) 

and a representative of its derivative T' by the embedding of the 5a distribution, 

1 / CL — r \ 

T'(r) = lim Um(5a) (r) = Umlim -r^f ). (9.3) 

a^O e^O a^O e^O e V £ / 

This implies that T has properties similar to Heaviside's step-function, and T' to those of 
Dirac's (5-function, with the fundamental difference that they are ^-functions which can be 
freely multiplied, differentiated, integrated, and combined with any C°° function. Moreover, 
the combination r^"T(r) can be given a precise meaning, and the number in Q, i.e., 
0(e'^), Vg G N, assigned to its point- value at r = 0. More precisely, the essential property of 
Tis 

1 1 f for r = 0, 

(_T)(r) = -T(r) = { (9.4) 
w > for r > 0, 
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in which the symbol r~"T(r) is interpreted as a whole such that the point- value r~"T(r)|^^Q = 
is well defined. This property has many practical consequencies: For example, it enables to 
write the potential 0(r) as in (19.11) . and it allows to freely associate the powers of r in products 
such as 

rP( -T{r)] f — T(r)) = -^T\r). (9.5) 

The detailed proof of (19.41) . which is given in [[T6ll . is not very complicated but somewhat 
lengthy. However, it is intuitively evident that the ^-identity (19.41) is simply the equivalent 
modulo O(e^), Vg G N, of the distributional identity 

1 1 f for r = 0, 

(-Hj(r) = -H„(r) = ^ ' (9.6) 

r" - r" - r " for r > 0, 

where the effect of is to 'cut-off' the singularity at r = as in (18.11 ). 

Equation (19.41) enables to work in Q with a much larger class of test functions than just the 
usual T(x ) = T(r, 9, 0) G ©(M^) of standard distribution theory. This class of 'generalized 
test functions' is 

jr .= ) = VT G V{R^), Vn G No}. (9.7) 

Due to (19.41) . the functions F(r, 9, 0) G JF(R^) are well defined at r = when multiplied by 
T or any of its derivatives. This enables to integrate the ^-functions T/r", T'/r", etc., which 
arise when calculating the electromagnetic fields of point charges and their derivatives, as well 
as any algebraic combination of them lfT6l . We will however content ourselves with just a few 
basic integration formulas, which are proved in the Appendix, i.e., 

) /"OO 

drT{r)F{r) = lim / dr F{r), (9.8) 







oo 

2/ 







drT\r)F{r) = lim / dr F{r), (9.9) 



oo 



drV{r)F{r) =limF(a), (9.10) 







1 

dr T{r)T\r)F{r) = lim -F(a), (9.1 1) 

a— >0 2 

r dr{rf{r)T{r) = limqo]^ + qi]T'(0), (9.12) 

where F{r) G !F, and T(r) G P. As expected, the ^-functions T and T' have properties 
similar to the distributions H and 5. Also, whereas formulas (19.81 - 19.1 11) do not depend of the 
particular representative of T (i.e., on the shape of r] provided r]{x) dx = 1), formula 
(19.121) explicitly depends on it because 

/ + 00 f + oo 

dxr]'^{-x), and Cii] = / dxxrj^{-x), (9.13) 
-oo J — oo 

where C[i] = if is even. This kind of undeterminedness is an intrinsic feature of the 
nonlinear context: Products of ^-functions lead to results which in general depend on the 
form of their representatives, a form which is determined by the physical problem. 
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To conclude this section we recalculate the Coulomb field and charge density in the 
T-formalism, that is starting from the potential (I9.1L which being a C°° expression permits to 
calculate as in elementary vector analysis. Thus, since Vr = u, 

E(f) = -V0(r) = e(^T(r) - -T(r))u, (9.14) 
Vr^ r / 

which is fully equivalent to (|8.7|) . The calculation of p(r) is also elementary, and leads to a 
greatly simplified result. Indeed, as V ■ m = 2/r, 

Anpir) = V ■ E{f) = e(lT(r) - ^-T'ir)) ^ 

+ e(-^T{r) + \r{r) + \r{r) - -T"(r)) = -e-T"(r),(9.15) 

which is much simpler than (18.101 ), and is easily seen to be associated to the usual three- 
dimensional charge-density because T"(r)/r x — (5(r)/r^ in R''. Moreover, this expression 
has the virtue of clearly showing the 'origin' of the charge density: The T(r) factor in the 
potential (I97TT ). 



10. Point-charge self-energy in V and Q 

Now that we have derived the Colombeau generalized functions corresponding to the Coulomb 
potential, field, and charge distributions, it is of interest to show that the 'extra' 5-like terms 
in (18.71) and (19.141) — which disappear when considering these distributions in V rather than 
in ^ — are physically significant. To do this we calculate the self-energy of a point charge 
using the electric field E defined according to three theories: The classical theory, distribution 
theory, and the Q theory, but using in all three cases the same self-energy expression derived 
from the Maxwell energy-momentum tensor. That is, in the classical theory in which E{f) is 
just the Coulomb field er/r^, the integral 

UseM ■■= III d^r E^ = l I dr r^^ = ^ lim - = oo. (10.1) 



8vr JJJ^3 2 Jo 2 r^o r 

In distribution theory we take for the Coulomb field the distribution E{f) = eRgf/r^ 
defined by (18.81) . Then, apart from expressing the self-energy Usc\i{T) as a function of the 
cut-off a, we still have the same divergent result 

t^seif(l) = = ^ /// d'rP = |lim - = oo, (10.2) 

Stt Svr JJJ^i 2 a-*o a 

even if E'^ is evaluated on a test-function T ^ 1. Thus, whereas the distribution (18.81) is 
meaningful for all r > 0, and gives sensible results for expressions linear in E evaluated on 
any test-function, it does not give a sensible result for the self-energy, which is quadratic in E. 
In particular, it is not possible to take the limit a — > which is mandatory for having a point 
charge. 

— * 

We now calculate the self-energy in Q, where the square of E^ is well defined. With the 
Coulomb field expressed as (19.141) . the self-energy is 

^7seif = ^ /// d'rE^ = i- r drr^(-J{r)--r{r}' 
Stt .///roa 2 . L r 
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= 2" /o l;^^'^'') ~ -T(r)r (r) + (r)^(r)) , (10.3) 

where all terms could freely be multiplied, and the factor simplified with the 1/r" factors, 
because everything is C°°. Then, after integration, the first two terms cancel each other exactly. 
Indeed, due to the identity 

-r\r)]' = i-T2(r) - -T(r)r(r), (10.4) 



we get, integrating (110.31) by parts. 



t4elf — ^ 

r 



2/ 



e 



+ - / rfr (r)'(r), (10.5) 



where the first term is zero on account of (|9.4|) . whereas the integral of the (T')^ term gives 
by means of (19.121 ) the result 

g2 I r+co 

t4eif = 77 lim - / dxrf{-x). (10.6) 

Because of the cancellation of the first two terms in (110.31) . only the square of the T' function, 
whose support is a point, contributes to the self-energy. This means that while the self-energy 
is infinite in the limit e ^ 0, this energy is now 'concentrated' at the location of the point 
charge rather than spread over the whole space surrounding it. 

Summarizing, when the self-energy is calculated in Q rather than in V , the divergent 
classical Coulomb-field self-energy (110.21) is canceled by the mixed term in the integral (110.31) . 
and that cancellation is exact and independent of the limiting sequence a ^ which is implicit 
in the symbol T. The sole contribution to the self-energy comes then from the (5^(r) term in 
that integral. This yields the result (110.61) which depends only on the shape of the moUifier 
7] and on the regularization parameter e, and which may be renormalized to a finite quantity 
such as the mass of the point charge. 

We have therefore obtained the physically remarkable result that in the Colombeau algebra 
— in which the multiplication of distributions is a meaningful operation — the self-energy of 
a point-charge is entirely located at the position of the charge, and solely due to the square 
of the T'(r) term in the electric field (|9.14l) . which itself derives form the T(r) factor in the 
potential (lOT ). 



11. Application of T-methodology to linear problems 

The final purpose of this paper is to confirm the methods used in papers [[m[T2| . in which only 
linear problems were considered. This is easily done by referring to the foregoing subsections, 
and by making the simplifications that are possible in that context. 

Indeed, in such applications there are no products of distributions such as T^, TT', or 
(T')^. The only remaining integration formulas are (19.81) and (19.101) . which reduce to the 
usual equations defining the properties of the Heaviside and Dirac distributions provided 
one substitutes T(r) H(r) and T'(r) — > 5{r). However, as was noted in [[TT|. the 
property T(r)|r=o = T'(O) = is necessary for the consistency of the formalism so that we 
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keep the notation T for that generalized function. On the other hand, there is no absolute 
necessity to distinguish between T'(r) and 5{r). Thus, instead of T"(r), T'"(r), etc., one can 
systematically use the associated distributional expressions ~5{r)/r, 25(r)/r^, etc., that is, 
(T')(") X 5(r)/r". 

Then, after all these simplifications, one may wonder why it is necessary to refer to 
generalized functions to justify the T-formalism in the linear context. The answer is that the 
simple rules introduced and used in IfTTl [T2l — which imply assigning the point value to 
T(0) and working with T and S as if they were C°° functions — only make fully sense in a 
framework of generalized functions such as a Colombeau algebra. For this reason, while it is 
possible to forget about most of the technicalities of that theory when using the T-formalism in 
linear problems, it is important to appreciate that working with T' instead of 5, as in Secs.[9]-fT0l 
of the present paper, is not much more complicated, and possibly less prone to mistakes. 

Finally, and to conclude, it is perhaps important to stress that the Colombeau formalism 
highlights the non-unicity of the 'microscopic' representations which at the end of a calculation 
give a physically meaningful result at the 'macroscopic' level, i.e., when evaluated on a 
test- function. This is why the methodologies that we labeled 'standard' and 'T' are most 
probably not the final theories which at the 'microscopic' level have a physically meaningful 
interpretation of their own: This is possibly much more the case of the methodology first 
introduced by Frank R. Tangherlinij^ which has the virtue of associating the ^-function T 
to the discontinuity of a truly fundamental quantity, the absolute value \ f\ of the distance 
between a source-point and a test-point, rewritten as \ f\ = rT(r) so that this discontinuity is 
properly taken into account [J7il . 
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13. Appendix: Proof of integration formulas 

To prove (19.81) and (|9.9I) it suffices to refer to (19.41) which specifies that apart from the point 
r = the function T(r) can be identified with one. Thus, for any m, 

oo POO 

dr T™(r) F(r) = Um / dr F{r). (13.1) 







To prove (19.101 ) we integrate by parts its left-hand side, i.e., 

/ dr r(r)F(r) = T(r)F(r) - / dr T(r)F'(r) (13.2) 
Jo Jq 

POO 

= -\im dr F'{r) = lim F (a), (13.3) 
where (l94l) and (fTO) were used in (fTT2l) . and F(oo) = in (fT33l) . 

^ This method has been independently rediscovered by the author and by others. See ifTTI and references therein. 
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Similarly, the identity (T^)' = 2TT' enables to integrate (19.1 II ) by parts, 



dr T(r)T'(r)F(r) = -T^{r)F{r] 



1 

2 Jo 



dr T^(r)F'(r) 



(13.4) 



lim- / dr F'{r) = \im-F{a), (13.5) 



which, using (|94l) and (|9^ in (|13.4D . proves (19.111) because T(oo) = 0. 

Finally, for the integration formulas of (T')^, i.e., (I9.12j - [97T31 ). we refer to Eq. (16.31) in 
which the relevant calculations are made. 
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